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2.  Linear  Systems  and  Optimal  Control,  (with  G.  Chen),  Springer  Series  in  Infor¬ 
mation  Sciences  #  18,  Springer- Verlag,  1988. 

3.  Multivariate  Splines,  CBMS-NSF  Series  in  Applied  Mathematics  #  54,  SIAM, 
Philadelphia,  1988. 

Papers 

1.  Image  reconstruction  by  bivariate  quadratic  splines,  (with  A.K.  Chan  and  K.B. 
Chan),  IEEE  Trans,  ASSP,  Vol.  36  (1988),  1525-1529. 

2.  A  natural  formulation  of  quasi-interpolation  by  multivariate  splines,  (with  H. 
Diamond),  Proc.  Amer.  Math.  Society.  Vol.  99  (1987),  643-646. 

3.  Interpolation  by  multivariate  splines,  (with  H.  Diamond  and  L.A.  Raphael), 
Math.  Comp.,  Vol.  51  (1988),  203-218. 

4.  On  multivariate  vertex  splines  and  applications,  (with  M.J.  Lai),  in  Topics  in 
Multivariate  Approximation ,  Academic  Press,  N.Y.,1987,  pp. 19-36. 

5.  Bivariate  C 1  quadratic  finite  elements  and  vertex  splines,  (with  T.X.  He),  Math. 
Comp,  Vol.  54,  No.  189  (1990),  169-187. 

6.  Shape-preserving  Quasi-interpolation  and  interpolation,  (with  H.  Diamond  and 
L.A.  Raphael),  J.  Comp,  and  Appl.  Math.,  Vol.  25  (1989)  169-198. 

7.  Shape-preserving  interpolation  by  bivariate  C 1  quadratic  splines,  (with  H.C.  Chui 
and  T.X.  He)  CAT  Report  #  148. 

8.  Mutlivariate  vertex  splines  and  finite  elements,  (with  M.J.  Lai),  J.  Approx.  The¬ 
ory,  Vol.  60,  No.  3  (1990),  245-343. 

9.  On  bivariate  super  vertex  splines,  (with  M.J.  Lai),  Constr.  Approx.  Vol.  6 
(1990),  399-419. 

10.  Modified  extended  Kalman  filtering  and  a  real-time  parallel  algorithm  for  system 
parameter  identification,  (with  G.  Chen  and  H.C.  Chui),  IEEE  Trans.  Auto  Cont, 
Vol.  35,  No.  1  (1990),  100-104. 

11.  System  reduction  via  truncated  Hankel  matrices,  (with  X.  Li  and  J.D.  Ward), 
Math,  of  Control,  Signals,  and  Systems.  To  appear. 

12.  On  the  dimension  of  bivariate  superspline  spaces,  (with  T.X.  He),  Math.  Comp. 
Vol.  53  (1989),  219-234. 

13.  On  the  convergence  of  s-numbers  of  compact  Hankel  operators,  (with  X.  Li  and 
J.D.  Ward),  CAT  Report  #184,  Circuits,  Systems,  and  Signal  Processing.  To 
appear. 

14.  Approximation  and  interpolation  formulas  for  real-time  applications,  (with  H. 
Diamond),  Trans,  of  the  7th  Army  Conference  on  Appl.  Math,  and  Computing, 
1989,  765-772. 

15.  Vertex  splines  and  their  applications  to  interpolation  of  discrete  data,  in  Compu¬ 
tation  of  Curves  and  Surfaces ,  Ed.  by  W.  Dahmen,  M.  Gasca,  arid  C.A.  Micchclli, 
Kluwer  Academic  Publishers,  1990,  137-181. 
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1C.  Rate  of  uniform  convc  rgence  of  rational  functions  corresponding  to  best  approx  - 
imants  of  truncated  Trankel  operators,  (with  X.  Li  and  J.  D.  Ward),  Math  of 
Control,  Signals,  ana  Systems.  To  appear. 

17.  Construction  and  applications  of  interpolation  formulas,  in  Multivariate  Approx¬ 
imation  and  Interpolation,  Ed.  by  W.  Haussmann  and  K.  Jetter,  ISNM  Series  in 
Math.,  Birkhauser  Verlag  Basel,  11-23. 

18.  Solution  of  the  four-block  problem  via  minimum-norm  interpolation,  (with  X. 
Li),  CAT  Report  #214. 

19.  A  characterization  of  multivariate  quasi-interpolation  formulas  and  its  applica¬ 
tions,  (with  H.  Diamond),  Num.  Math,  Vol.  57  (1990),  105-121. 

20.  A  general  framework  for  local  interpolation,  (with  H.  Diamond),  CAT  Report 
#190,  Texas  A&M  University,  Numer.  Math.  To  appear. 

21.  A  cardinal  spline  approach  to  wavelets,  (with  J.  Z.  Wang),  CAT  Report  #211. 

22.  On  compactly  supported  spline  wavelets  and  a  duality  principle,  (with  J.  Z. 
Wang),  CAT  Report  #213. 

23.  A  general  study  of  maximal  robust  stability  regions,  (with  H.N.  Mhaskar),  Cir¬ 
cuits,  Systems,  and  Signal  Processing.  To  appear. 

24.  Radial  basis  function  approach  to  interpolation  of  large  reflecting  surface,  (with 
A.  Iv.  Chan  and  L.  T.  Guan),in  SPIE  Vol.  1251  Curves  and  Surfaces  in  Computer 
Vision  and  Graphics,  (L.  A.  Ferrari  and  R.  J.  P.  de  Figueiredo,  Eds.),  1990,  62-72. 

25.  Application  of  generalized  vertex  splines  to  boundary  element  method  for  elec¬ 
tromagnetic  scattering,  (with  A.  K.  Chan  and  T.  X.  He),  in  6th  Annual  Review 
of  Progress  in  Applied  Computational  Electromagnetic  Society,  1990,  329-337. 

26.  On  multivariate  robust  stability,  (with  H.  N.  Mhaskar),  CAT  Report  #216. 

27.  Curve  design  and  analysis  using  splines  and  wavelets,  Proc.  8th  Army  Conf.  on 
Appl.  Math,  and  Computing,  June  1990.  To  appear. 
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Q.  Liu 
J.  Wang 

Computer  programmers  (Undergraduate  students) 
Scott  Bowers 
Susan  Fojtasek 
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III.  Brief  outline  of  research  findings 

The  objective  of  this  report  is  to  give  a  very  brief  summary  of  the  research  findings 
on  the  project  under  Contract  Number  DAAL  03-87-K-0025  sponsored  by  SDIO/IST  and 
managed  by  the  U.S.  Army  Research  Office.  Technical  details  are  not  included  since  most 
of  them  have  been  reported  as  semi-annual  progress  reports.  This  project  covered  the 
period  February  1,  1987  -  July  30,  1990  which  includes  a  six-month  extension  at  no  cost. 
The  list  of  publicatins  in  Section  I,  page  1  -  page  2,  will  also  be  used  as  references.  Other 
references  are  listed  at  the  end  of  this  summary,  We  divide  our  results  into  five  categories 
as  follows: 

1.  Super  splines  and  vertex  splines. 

Our  results  in  this  area  are  found  in  [4, 5, 8, 9]  with  applications  to  interpolation  to  be 
discussed  in  the  next  section  and  other  applications  in  Section  4.  It  is  well  known  that 
every  spline  function  in  one  variable  can  be  represented  as  a  5-spline  series.  However, 
with  only  a  few  exceptions,  this  nice  local  representation  does  not  generalize  to  bivariate, 
and  more  generally  multivariate,  splines.  One  of  the  reasons  is  that  the  CT  piecewise 
polynomial  functions  on  a  triangulation  have  very  complicated  local  structures,  and  the 
one  that  is  crucial  to  the  computational  aspects  is  the  extra  free  parameters  on  lower- 
dimensional  manifolds  that  constitute  the  triangulation.  This  phenomenon  is  probably 
well  known  to  some  of  the  finite  element  people.  In  our  work  [4],  we  introduced  the  notion 
of  super  splines  in  order  to  drop  the  “useless”  parameters.  For  instance,  a  C 1  quintic 
bivariate  piecewise  polynomial  is  restricted  to  be  C 2  at  the  vertices.  In  our  paper  [9], 
we  considered  the  smallest  subspace  of  super  splines  in  S-ir+2  (i-e-  Cr  piecewise  bivariate 
polynomials  with  total  degree  3r  +  2)  and  showed  that  the  approximation  order  is  full, 
that  is  3r  +  3  as  shown  in  de  Boor  -  Hollig  [Si].  In  fact,  this  super  spline  subspace  has 
a  basis  consisting  of  vertex  splines;  and  hence,  a  spline  series  representation  is  obtained. 
The  notion  of  vertex  splines  was  introduced  by  Chui  and  Lai  in  [S9] .  Later  development  of 
super  and  vertex  splines  can  be  found  in  [5,8,9]  and  recent  work  of  Alfeld,  Hollig,  Ibrahim, 
Piper,  Schumaker,  and  others  (cf.  the  articles  in  [S10,S12,S13,S28]  and  references  therein). 

2.  Quasi-interpolation  and  interpolation  formulas. 

The  central  theme  in  this  area  of  research  is  construction  of  “real-time”  interpolation 
formulas,  and  real-time  implementation  of  such  formulas.  In  our  work  [2,3],  we  developed 
local  quasi-interpolation  schemes  by  introducing  the  so-called  “Neumann-series”  approach. 
Here,  quasi-interpolation  means  optimal-order  local  linear  approximation.  In  [19],  we  gave 
a  complete  characterization  of  all  quasi-interpolation  fornmulas.  The  most  important  ap¬ 
plication  in  [19]  is  construction  of  quasi-interpolation  using  arbitrary  sample  points  and 
various  data  information.  To  change  quasi-interpolation  to  actual  local  interpolation,  pre¬ 
serving  the  optimal  order  of  approximation,  we  introduced  in  [14,20]  a  (noncommutative) 
“blending”  method  by  using  vertex  splines  or  generalized  vertex  splines.  The  emphasis 
is  localncss  and  optimality  of  the  order  of  approximation.  Hence,  “table-look-up”  can  be 
used  for  real-time  interpolation;  or  equivalently,  an  FIR  (moving-average)  scheme  can  be 
implemented.  The  general  theory  is  extended  to  nonuniform  grid  and  an  application  to 
solving  nonlinear  PDE’s  is  discussed  in  [17],  and  many  other  applications  and  examples 
are  included  in  my  tutorial  article  [15],  written  for  the  NATO  Graduate  Studies  on  “Com¬ 
putations  of  Curves  and  Surfaces”.  In  our  work  [6,7],  we  studied  interpolation  formulas 
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that  preserve  the  shapes  of  the  data.  In  particular,  an  almost  real-time  method  is  given  in 
[6],  and  a  comprehensive  studr  m  shape-criteria  is  given  in  [7]. 

3.  General  multivariate  spline  theory. 

Our  articles  in  this  area  include  [4,5,7,8,9,12,19,20].  As  discussed  in  the  last  two  sec¬ 
tions,  the  notions  of  super  splines  and  vertex  splines,  and  the  methods  of  Neumann 
series  approach  and  noncommutative  blending  were  introduced  in  this  series  of  pa¬ 
pers.  To  facilitate  the  construction  schemes,  the  dimensions  of  the  super  spline  subspaces 
should  be  determined.  The  first  results  in  this  area  are  given  in  our  paper  [12],  where  a 
dimension  criterion  is  also  introduced.  The  general  theory  of  super  spline,  vertex  splines, 
quasi-interpolation,  and  interpolation  in  higher  dimensions  is  included  in  my  CBMS-NSF 
monograph  [S3].  In  the  second  edition  of  my  encyclopedia  article  [S2],  cardinal  spline 
interpolation  is  related  to  wavelets,  a  topic  to  be  discussed  in  Section  5. 

4.  Applications  to  engineering  problems. 

Our  results  in  this  area  are  included  in  [1,10,11,13,16,18,23,24,25,26].  In  addition  to 
construction  and  implementation  of  real-time  approximation  and  interpolation  schemes, 
we  have  found  almost  immediate  applications  to  image  reconstruction  [1],  design  of  large 
reflector  surfaces  [24],  and  an  application  to  the  boundary  element  method  for  electro¬ 
magnetic  scattering  [25].  On  the  other  hand,  real-time  methods  should  be  implemented 
with  the  Kalman  filter.  In  our  work  [10],  a  parallel  algorithm  is  developed  to  improve 
the  extedned  Kalman  filter  significantly  and  applications  to  real-time  system  parameters 
identification  are  given.  In  the  second  edition  of  our  monograph  [S7],  which  is  under 
preparation,  this  method  is  discussed  in  some  details.  When  frequency-domain  methods 
are  considered,  the  problem  of  stability  must  be  considered.  By  using  the  Hankel-norm 
approximation,  stability  is  guaranteed.  In  our  work  [11,13,16],  computational  efficiency 
and  accuracy  are  discussed  relative  to  truncation  of  the  Hankel  operator  and  the  estimates 
are  given  in  terms  of  s-numbers  (or  singular  values).  A  related  application  is  the  solution 
of  the  so-called  four-block  problem  in  systems  theory.  This  is  studied  in  [18].  Robust 
stability  are  studied  in  [23]  for  the  one-variable  (i.e.  SISO)  setting,  and  in  [26]  for  the 
multi-vanable  (i.e.  MISO)  setting. 

5.  Wavelets. 

Perhaps  the  most  exciting  development  in  this  project  is  the  introduction  of  cardinal- 
spline  wavelets  in  our  work  [21,22].  This  was  done  toward  the  end  of  the  extended  period  of 
this  project  which  overlaps  with  the  new  project  (under  Contract  No.  DAAL  03-90-0091). 
This  new  project,  with  the  same  title,  should  be  considered  as  continuation  of  the  old 
one  for  which  this  report  is  prepared.  The  multiresolution  analysis  introduced  by  Meyer 
[S22]  and  Mallat  [S19,  S20,  S21]  is  followed.  As  is  well-known,  Daubechies’  wavelets  [S16, 
S 1 7]  are  the  only  available  compactly  supported  orthonormal  wavelets.  However,  although 
there  are  efficient  algorithms  (cf.  [S 17] ) ,  these  wavelets  do  not  have  explicit  formulations 
and  are  certainly  not  suitable  for  real-time  applications.  In  addition,  since  they  are  not 
symmetric  or  anti-symmetric,  their  non-linear  phase  property  creates  distortion.  In  [21], 
we  introduced  the  notion  of  cardinal-spline  wavelets.  These  wavelets  are  constructed  by 
using  cardinal  fundamental  splines  of  Schoenberg.  Since  they  arc  explicitly  formulated, 
symmetric  for  even-order  splines,  and  antisymmetric'  for  odd-order  splines,  the,  are  readily 
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available  for  interpolation  and  wavelet  decomposition  of  data  in  real-time.  An  application 
to  curve-design  is  given  in  [27],  and  a  comprehensive  paper  on  application  to  real-time 
signal  analysis  is  under  preparation.  The  general  theory  of  our  new  wavelets  is  studied  in 
[S14],  partially  supported  by  the  new  project.  This  theory  includes  both  the  Daubechios 
wavelets  and  our  cardinal-spline  wavelets  as  special  cases.  A  fairly  complete  analysis  of 
the  minimally  supported  cardinal  spline  wavelets  (introduced  in  [22])  is  given  in  [S15]. 
A  survey  article  [S3]  has  been  written  and  will  be  published  by  Academic  Press,  and  a 
tutorial  article  [S4]  which  is  a  compilation  of  the  five-hour  lecture  series  I  gave  in  Lancaster. 
England,  last  July,  will  be  published  by  the  Oxford  University  Press. 
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